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ABSTRACT 

We compare two methods of distance determination to spiral galaxies using 
optical/near-infrared (NIR) observations, the (I-K) versus Mk colour - absolute mag- 
nitude (CM) relation and the /and if-band TuUy-Fisher relation (TFR). 
Dust-free colours and NIR absolute magnitudes greatly enhance the usefulness of the 
NIR CM relation as a distance indicator for moderately to highly inclined spiral galax- 
ies in the field (inclinations between ~ 80° and 90°); by avoiding contamination by 
dust the scatter in the CM relation is significantly reduced, compared to similar galaxy 
samples published previously. The CM relation can be used to determine distances to 
field spiral galaxies with Mk > —25.5, to at least Mk ~ —20. 

Our results, supplemented with previously published observations for which we can 
- to some degree - control the effects of extinction, are consistent with a universal 
nature of the CM relation for field spiral galaxies. 

High-resolution observations done with the Hubble Space Telescope can provide a 
powerful tool to calibrate the relation and extend the useful distance range by more 
than a factor of 2 compared to ground-based observations. 

The intrinsic scatter in the NIR CM relation in the absolute /f-band magnitudes is 
~ 0.5 mag, yielding a lower limit to the accuracy of distance determinations on the 
order of 25%. 

Although we find an unusually low scatter in the TFR (probably a statistical accident), 
a typical scatter in the TFR would yield distances to our sample galaxies with uncer- 
tainties of only ~ 15%. However, one of the main advantages of the use of the NIR 
CM relation is that we only need photometric data to obtain distance estimates; use of 
the TFR requires additional kinematic data, although it can be used to significantly 
greater distances. 

Key words: distance scale - galaxies: photometry - galaxies: spiral ~ galaxies: statis- 
tics - infrared: galaxies 



1 PHOTOMETRIC DISTANCE INDICATORS 

The study of highly- inclined disc-dominated galaxies offers 
the unique opportunity to avoid the confusing effects of in- 
terstellar extinction on the optical and near-infrared (NIR) 
appearance of their dominant (old) stellar disc population. 
In this paper we use NIR absolute magnitudes and I-K 
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colours that were determined in those regions of our sta- 
tistically complete sample of edge-on field galaxies that are 
least affected by extinction, to derive a "dust-free" colour 
- absolute magnitude (CM) relation. Dust-free colours and 
NIR absolute magnitudes greatly enhance the usefulness of 
the NIR CM relation as a distance indicator for moderately 
to highly inclined spiral galaxies in the field. 

Alternatively, large galaxy surveys in the NIR facili- 
tate the use of the NIR TuUy-Fisher relation (TFR; TuUy 
& Fisher 1977) as an accurate tool to obtain distances to 
spiral galaxies in clusters. In this paper we discuss the / and 
_R'-band TFRs derived from our sample galaxies and the ac- 
curacy with which they can be used to estimate distances to 
field spiral galaxies. 
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2 THE NIR COLOUR-MAGNITUDE 
RELATION 

The tightness of the CM relation for early-type galaxies (as 
first established by Baum [1959] and de Vaucouleurs [1961]), 
makes it potentially useful as a distance indicator, as was 
first suggested by Sandage (1972). In this paper we inves- 
tigate whether the NIR CM relation is also useful as a di- 
agnostic tool to estimate distances to spiral galaxies in the 
field. 



2.1 Elliptical versus spiral galaxies 

Visvanathan & Griersmith (1977) extended the range of 
galaxy types from elliptical/SO galaxies to early-type spi- 
rals (so/a to Sab), and found (within the errors) exactly 
the same optical CM relation for the early-type Virgo clus- 
ter spirals as had been found for E/SO galaxies, but with a 
larger excess scatter (see also Visvanathan & Sandage 1977; 
Griersmith 1980; Bower, Lucey & Ellis 1992b; Peletier & de 
Grijs 1998). 

Despite the good agreement between the elliptical and 
spiral galaxy CM relations, later-type spiral galaxies occupy 
a different region in the CM diagram (e.g., TuUy, Mould & 
Aaronson 1982; Mobasher, EUis & Sharpies 1986; Peletier 
& de Grijs 1998). Apart from this distinction between the 
loci of late-type spirals and E/SO's, Griersmith (1980) also 
noticed that differences in zero point of the CM relations 
seem to follow a systematic trend along the Hubble sequence: 
colours become systematically bluer for later Hubble types 
(see also Peletier & de Grijs 1998), although the slopes for 
the individual galaxy types are the same (within the errors) 
compared to each other and to E/SO galaxies. 



2.2 A universal CM relation? 

Following the discovery of the CM effect for the early-type 
Virgo cluster members, many attempts have been made to 
unambiguously determine the universality of the relation- 
ship (e.g., Visvanathan & Sandage 1977; Visvanathan & Gri- 
ersmith 1977; Griersmith 1980; Aaronson, Persson & Frogel 
1981; Bower et al. 1992a,b). Provided that the CM relation 
is universally applicable it can be used as a distance indi- 
cator. In general, the UV-optical CM relation for early-type 
galaxies was found to be universal and independent of envi- 
ronment to a high degree (for a recent review see Ellis et al. 
1997). However, the tests done to study the universality of 
the CM effect were challenged by a number of studies (e.g., 
Faber 1977; Burstein 1977; Larson, Tinsley & CaldweU 1980; 
Aaronson et al. 1981, and more recently Abraham et al. 1996 
and van Dokkum et al. 1998), claiming that non-negligible 
environmental effects were playing a role in the observational 
evidence. To undertake an independent study of the univer- 
sality of the CM relation. Bower et al. (1992a,b) obtained 
new observations of Virgo and Coma cluster galaxies. Their 
observations support the argument that a universal CM re- 
lation for cluster E/SO galaxies, over the entire wavelength 
range, is likely, although the dispersion is expected to be 
considerable. Bower et al.'s (1992b) photometry allows the 
CM relation for early-type cluster galaxies to be used to 
estimate distances accurate to ~20% per galaxy. 



In this paper we will use our own optical and NIR obser- 
vations of edge-on spiral galaxies, supplemented with similar 
samples from the literature, to investigate the possible uni- 
versality of the CM relation for spiral galaxies in the field. 



2.3 Advantages of NIR observations 

For the E/SO galaxies in the Virgo cluster the change of 
colour with absolute magnitude is greatest in the ultraviolet, 
and decreases significantly towards redder wavelengths (Vis- 
vanathan & Sandage 1977, and references therein). However, 
the CM effect shows up again in the optical-NIR regime, in 
the sense that V-K is bluer for intrinsically fainter galaxies 
(e.g., Aaronson et al. 1981, TuUy et al. 1982). The usefulness 
of NIR observations of spiral galaxies in the field for measur- 
ing extragalactic distances depends on two main questions 
(e.g., Aaronson et al. 1981): 

(i) Do spiral galaxies in the field follow a universal CM 
relation in the NIR; 

(ii) If so, is the scatter sufficiently small for a useful ap- 
plication of the relation? 

The main advantages of NIR observations compared to 
observations in the blue passbands are their relative insen- 
sitivity to contamination by the presence of young stellar 
populations and dust. The absorption corrections for dust in 
external galaxies, which are largest in the blue, are difficult 
and controversial; therefore, we were motivated to study the 
CM relation for our sample of edge-on spiral galaxies in the 
NIR, using I-K colours and absolute _R'-band magnitudes. 

In Peletier & de Grijs (1998) we find that the scatter in 
the NIR CM relation for (field) spiral galaxies can entirely be 
explained by observational uncertainties. Moreover, we show 
that the slope of the relation is steeper for spirals than for 
ellipticals. We can explain this if we assume that the CM 
relations for spiral and elliptical galaxies are intrinsically 
different, in the sense that the stars in spiral galaxies are 
younger than those in ellipticals and the fraction of young 
stars in a spiral galaxy (i.e., its "age") is determined solely 
by the galaxy's luminosity, and not by its environment. The 
CM relation formed by elliptical galaxies, on the other hand, 
is generally attributed to changes in metallicity (Peletier & 
de Grijs 1998, and references therein). 

In Sect. M we discuss the nature of our observations and 
the corrections that have to be applied for an accurate study 
of the NIR CM relation. We supplement our data with re- 
cently published / and A'-band data in Sect. 0. In Sect. 5 
we discuss the implications of our results and look into the 
applicability of a tight NIR CM relation using HST obser- 
vations. We discuss the accuracy of the NIR TFR in Sect. 
6, and we give a summary of our results and conclusions in 
Sect. 7. 



3 CONSTRUCTING THE CM RELATION 

The observations of and the reduction techniques applied to 
our sample of highly-inclined spiral galaxies, as well as the 
resulting photometric accuracy, have been described in de 
Grijs (1997, 1998). Here, we summarize our selection crite- 
ria, applied to the galaxies in the Surface Photometry Cat- 
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alogue of the ESO-Uppsala Galaxies (ESO-LV; Lauberts & 
Valentijn 1989): 

• Inclination selection: i > 87°; 

• Diameter selection: D^ > 2. '2 (-D;^: blue angular size 
at a surface brightness Es ~ 25 mag arcsec" ); 

• Galaxy type selection: SO - Sd, and 

• Morphology selection: they should be non-interacting 
and undisturbed. 

Although our NIR observations are relatively insensitive 
to disturbing dust effects (see, e.g., de Grijs, Peletier & van 
der Kruit 1997) , we need to apply several corrections to the 
observed quantities before they can be used to study the 
CM relation. 



3.1 Colour gradients and Galactic foreground 
extinction 

Colour gradients in edge-on galaxies have not been studied 
extensively. In a few well-studied edge-on spiral galaxies, 
colour gradients parallel to the major axis were found to be 
negligible (e.g., Hamabe et al. 1980, Jensen & Thuan 1982, 
van der Kruit & Searle 1982), or show an increasingly blue 
disc population with increasing radius (e.g., Sasaki 1987, 
Aoki et al. 1991), although such gradients are generally 
small. 

Vertical colour gradients in the discs of our sample edge- 
on galaxies are generally small or negligible (de Grijs et al. 
1997, de Grijs & Peletier, in preparation); in fact, the ver- 
tical colour gradients are smaller than or of the same order 
as the observational errors for our I-K colour determina- 
tions. The only statistically significant vertical colour varia- 
tion is the reddening due to dust in the galactic planes. As 
we showed in de Grijs & van der Kruit (1996) and de Grijs & 
Peletier (1997), the vertical scale height does not vary signif- 
icantly as a function of wavelength, which indirectly shows 
the absence of any significant vertical colour gradient. 

In de Grijs et al. (1997) we show that, once away from 
the dust lane, the vertical I-K colour profiles are generally 
(approximately) fiat and featureless, indicating that the ex- 
cess extinction in the / band (compared to that in the K 
band) is negligible at these heights above the galactic planes. 
The colours we use in this paper were determined from this 
(approximately) fiat part of the vertical colour profile at the 
minor axis (at 1.5hz < \z\ < 3.5/iz, where h^ is the expo- 
nential scale height), to avoid the reddening caused by the 
in-plane dust. Since young stellar populations are generally 
confined to regions close to the galactic planes, an additional 
advantage of determining colours away from the planes is 
that the contamination by emission from these young stel- 
lar populations is thus greatly reduced or negligible. 

Peletier & Balcells (1997) showed, for their sample of 30 
field spiral and lenticular galaxies, that the bulge colours on 
the minor axis and the inner disc colours taken in wedge 
apertures at 15° from the major axis at 2 Tf-band scale 
lengths are very similar. Therefore, the colour of an edge-on 
disc-dominated galaxy determined in a dust-free region at 
the minor axis can be considered as representative for the 
galaxy's dominant (old) stellar population, under the as- 
sumption that vertical colour gradients in the old-disc pop- 
ulation are smaller than the observational uncertainties. 

To deal with Galactic reddening, we used the Galactic 



extinction values in the / and K' bands, Ag,i and Aqj^i 
given by Schlegel, Finkbeiner & Davis (1998), thereby as- 
suming that the Galactic extinction can be approximated by 
a foreground dust screen. These full-sky dust maps are twice 
as accurate as the older reddening estimated by Burstein & 
Heiles (1978, 1984) in regions of low and moderate redden- 
ing and likely significantly more accurate in regions of high 
reddening (see also Hudson 1999). The Galactic extinction 
in the / and K' bandpasses was obtained assuming a stan- 
dard Rv = 3.1 extinction law (Schlegel et al. 1998, their 
Appendix B). 



3.2 Cosmological corrections and distance 
calibration 

The corrections for the effects of redshift (Tf- corrections) are 
generally small in the NIR. For our sample galaxies the K- 
corrections for the I-K colours range from 0.00 - 0.04 mag 
(e.g., Schneider, Gunn & Hoessel 1983), depending on galaxy 
type (e.g., Coleman, Wu & Weedman 1980), see Table |l|. 

Heliocentric velocities for the majority of our sample 
galaxies were obtained by Mathewson, Ford & Buchhorn 
(1992) and Mathewson & Ford (1996), which provides us 
with a homogeneous data set to base our absolute magnitude 
calculations on. To obtain absolute magnitudes we applied 
the formula for the systemic velocities adjusted for the solar 
motion with respect to the centroid of the Local Group given 
by Richter, Tammann & Huchtmeier (1987, see also Schmidt 
& Boiler 1992): 



VLG = «0 -I- Aw , 

where 

A« = -49.59 cos £ cos b -I- 306.95 sin £ cos b 
- 18.59 sin b 



(1) 



(2) 



In Table hi we give an overview of the observed and 
derived quantities used for the study of both the NIR CM 
and TF relations. 



4 COMPARISON WITH LITERATURE 

SAMPLES 

In Fig. hla we have plotted the CM relation derived from our 
"dust-free" I-K colours (determined in the region 1.5hz < 
zl < 3.5ftz) and the absolute .ftT-band magnitudes, M%, both 
corrected for Galactic foreground extinction. To assess the 
importance of contamination by dust, in Fig. hib we compare 
the NIR CM relation resulting from the use of integrated 
galaxy colours to that derived from the dust-free colours. It 
is clear that the importance of dust should not be under- 
estimated: the integrated colours are systematically redder 
than the "dust-free" colours. 

Careful examination of Fig. nk reveals a rather pecu- 
liar distribution of the absolute ii'-band magnitudes of our 
sample galaxies, into two clumps and a single, bright galaxy. 
We believe that this is not due to any selection effects other 
than the availability of observing time with NIR arrays. The 
corresponding distribution of absolute /-band magnitudes, 
which were obtained for about twice as many of our sam- 
ple galaxies, does not show any significant subclumping, al- 
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Table 1. Basic properties of the sample galaxies 

Columns: (1) Galaxy name (ESO-LV); (2) heliocentric velocity (Mathewson et al. 1992); (3) velocity correction for motion w.r.t. the 
centroid of the Local Group (Richter et al. 1987); (4)— (5) and (6)-(7) apparent /and ii'-band magnitudes and observational errors; (8) 
Galactic extinction in K'; (9) I-K colouv excess (Schlegel et al. 1998); (10) and (11) I-K color, corrected for Galactic foregound extinction, 
and observational error; (12) if-correction (interpolated from Schneider et al. 1980); (13) absolute iC-band magnitude (Hq = 100 km s~^ 
Mpc~^), corrected for Galactic foreground extinction and the effects of redshift; (14) peculiar velocity with respect to the centroid of 
the Local Group (Mathewson et al. 1992). 
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though the two brightest galaxies are ^ 1 mag brighter than 
the third brightest galaxy observed as part of our project. 

In Fig. y we compare our results with those of other 
galaxy samples published recently. Only very few samples 
with high-quality if-band observations have been published 
to date: 

(a) The sample of Bershady et al. (1994) and Bershady 
(1995) consists of 171 field galaxies of all types, selected 
from three high Galactic latitude fields (where the Galactic 
extinction is low), of which a subset of 143 is statistically 
representative in its sampling of the apparent colour distri- 
bution of galaxies. The median redshift of the sample galax- 
ies is z = 0.14. Absolute magnitudes were determined using 
luminosity distances, corrected for cosmological effects. 

(h) De Jong's (1996b) sample is a diameter- limited sample 
of 86 (nearly) face-on spiral galaxies. His absolute magnitude 
calculations are based on heliocentric velocities, corrected 
for Virgocentric infall using the model of Kraan-Korteweg 
(1986; de Jong, priv. comm.), and taking into account 
the Galactic foreground extinction estimates of Burstein & 
Heiles (1978, 1984). 

(c) The data presented by Andredakis, Peletier & Balcells 
(1995) and Peletier & Balcells (1997) consists of 37 field disc 
galaxies of types SO to Sbc, uniform in orientation on the 
sky, for 30 of which NIR observations are available. To ob- 
tain absolute magnitudes, they used the Galactic standard- 
of-rest corrections given by de Vaucouleurs et al. (1991). 



Corrections for foreground extinction were only applied to 
their /-band observations, since the K-hanA corrections were 
smaller than their observational errors. 

(d) TuUy et al.'s (1996) sample was taken from the Ursa 
Major cluster; they present / and /f'-band observations of a 
magnitude-limited sample of 70 disc-dominated galaxies. 

After scaling these and our observations to the same 
cosmology (using Ho = 100 km s~^ Mpc~^), we find that 
the most significant deviations from our CM relation are 
exhibited by samples (a) (Bershady et al. 1994, Bershady 
1995), and (d) (TuUy et al. 1996). 

In either of these samples the reddening due to dust 
is unknown. Bershady et al. (1994) and Bershady (1995) 
corrected their integrated / and /f-band colours neither for 
Galactic nor for internal extinction in their sample galaxies. 
Any comparison between our dust-free and their integrated 
colours will therefore show an additional discrepancy due 
to extinction effects. Moreover, their redshift distribution is 
different from those of the other samples, so that we also 
have to take into account possible evolutionary effects when 
comparing these relations. 

TuUy et al. (1996) corrected their photometry for both 
Galactic foreground and internal extinction, using TuUy 
& Fouque's (1985) inclination corrections. Although this 
method is to first order useful for low and intermediate incli- 
nations, there are severe problems for highly-inclined galax- 
ies. Moreover, for the analysis of a combined sample of spiral 
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Figure 1. (a) Correlation between absolute _ft'-band magnitude, M^, and dust-free I-K colours. The best-fitting NIR CM relation is 
indicated by the dashed line, for which the parameters are given in Table H (b) Corresponding CM relation derived from the integrated 
galaxy colours. Almost all data points are located redwards from the best- fitting CM relation shown in Fig. hla, due to internal interstellar 
extinction. All observational quantities have been corrected for Galactic foreground extinction. 



galaxies, extinction effects should be treated uniformly. In 
addition, Tully et al.'s (1996) sample was drawn from a clus- 
ter population, whereas the other samples consist (mainly) 
of field galaxies. Until it has been established unambiguously 
that environmental effects do not play a role in the deter- 
mination of the CM relation for spiral galaxies, we cannot 
use a mixture of field and cluster galaxy samples to base 
our conclusions on. Therefore, we will consider field galaxy 
samples only in our analysis of the NIR CM relation in this 
paper. 

Sample (b) (de Jong 1996b) consists of disc-dominated 
(nearly) face-on galaxies. The face-on orientation ensures 
that the effects of internal extinction on the integrated I- 
K colours are small (but not necessarily negligible) , as op- 
posed to more inclined galaxies, since less light is obscured 
by the dust lane, because the light dominating these colours 
is emitted by the old-disc stellar population located in front 
of the dust component (which is concentrated in the galactic 
plane region). De Jong (1996a) argues that the colours and 
colour gradients observed in his sample of face-on galaxies 
are determined by intrinsic physical processes in the discs 
rather than by extinction effects. However, it is likely that 
the larger scatter in his CM relation compared to ours is at 
least partly due to extinction effects affecting his integrated 
colours (which are likely greater for these colours than for 
our "dust-free' colours, due to line-of-sight integration ef- 
fects), in addition to the effects of colour gradients across 
his face-on discs, which vary as a function of galaxy type. 

Finally, sample (c) (Andredakis et al. 1995, Peletier & 



Balcells 1997) consists of early- type disc galaxies (SO-Sbc), 
of which they studied the disc and bulge components sepa- 
rately. From the orientation of the individual galaxies they 
obtained dust-free colours on the non-obscured side of the 
galactic planes. In addition, Peletier & Balcells (1997) as- 
sessed the importance of internal extinction by studying the 
change in the position of the galactic centres as a function of 
wavelength. Such a shift occurs when a dust lane in front of 
the true centre obscures more light on one side of the centre 
than on the other, causing the observed centre to shift as 
a function of optical depth, or passband. Their centre shift 
is in general within the typical errors, which provides addi- 
tional support to the assumption that extinction does not 
play a major role in their colour determinations. 

In Table H we compare the parameters of the NIR CM 
relation derived from our sample of 22 edge-on disc galaxies 
with both the CM parameters derived from the comparison 
samples in Fig. bl and those obtained from supplementing 
our data with samples (b) and (c), the "composite" CM 
relation, for which we are fairly confident that extinction 
effects are small. 

One might wonder to what extent the slope and zero 
point of our CM relation are determined by the brightest, 
outlying galaxy. We redetermined these parameters for our 
sample, excluding this brightest member, and found essen- 
tially identical results (see Table |2[ for the sample "without 
brightest"). This was of course to be expected, in view of the 
good agreement between our data on the one hand, and the 
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slope = -0.243 (±0.012; 
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(b) 
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Figure 2. Comparison of our Mx vs. (I-K) CM relation with previously published data. The filled circles show our sample galaxies; 
the open circles represent previously published data. The best-fitting NIR CM relations for the literature samples are indicated by the 
solid lines; for comparison, the dashed lines represent the mean CM relation derived from the data presented in this paper, (a) Bershady 
et al. (1994), Bershady (1995); (b) de Jong (1996b); (c) Andredakis et al. (1995), Peletier & Balcells (1997); (d) TuUy et al. (1996). 
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Table 2. Comparison of the NIR CM parameters: 

(I -K) = ahlK + b 

Columns: (1) Sample; (2) and (3) CM slope, a, and its error; (4) 

and (5) Intercept, 6, and its error; (6) Scatter in the CM relation. 



on dust-free colours may in principle be useful as a diagnos- 
tic to estimate distances. It provides independent distance 
estimates based on observational parameters. 



Sample 




a 


± 


b 


± 


0- (mag) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(a) Bershady 




-0.243 


0.013 


-3.48 


0.30 


0.256 


(b) de Jong 




-0.226 


0.024 


-3.38 


0.55 


0.222 


(c) Peletier 




-0.226 


0.042 


-3.63 


1.04 


0.173 


(d) TuUy (all) 




-0.193 


0.018 


-2.46 


0.37 


0.242 


(d) TuUy (T > 


4; 


-0.199 


0.025 


-2.67 


0.51 


0.256 


our 




-0.156 


0.021 


-1.99 


0.48 


0.133 


without brightest 


-0.157 


0.024 


-1.99 


0.55 


0.137 


composite 




-0.189 


0.013 


-2.59 


0.31 


0.205 



similar samples used to construct the composite CM relation 
on the other. 

The dashed line in Fig. H shows the composite NIR CM 
relation. The best-fitting CM relation derived from our data 
alone is shown by the dotted line. The close agreement be- 
tween both CM relations supports our assumption that we 
can indeed compare our observations to samples (b) and (c). 
In fact, the good agreement between our sample and samples 
(b) and (c) leads us to note that the NIR CM relation for 
the old-disc population of spiral galaxies in the field appears 
to be universal (see Peletier & de Grijs 1998). 

In Table we also give the standard deviations of the 
observed data points with respect to the best-fitting linear 
CM relation, which are indicative of the scatter in these 
CM relations. The best-fitting relations were determined by 
means of a bivariate minimization algorithm. It is immedi- 
ately clear that the scatter in the CM relation determined 
from our dust-free data points is of order 40% smaller than 
for both samples (b) and (c). 

Only if we can put firm constraints on the scatter in 
the CM relation for field spiral galaxies, it may be useful as 
a distance indicator. Therefore, we compared the scatter in 
the CM relation that we obtained from our dust-free disc 
colours to that from previous studies. 

From a detailed photometric study of Virgo and Coma 
cluster members. Bower et al. (1992b) found that for ellip- 
tical galaxies the scatter in the CM relation is dominated 
by the observational errors. The inclusion of SO and early- 
type spiral galaxies into their sample increases the observed 
scatter in the optical-NIR ( V-K ) CM relation. To be able 
to compare our results to those of Bower et al. (1992b), we 
applied the least-squares fitting technique we used to their 
data, and found close matches between our and their values 
for the r.m.s. scatter. Therefore, we can directly compare the 
scatter estimates obtained from the samples listed in Table 
||to Bower et al.'s (1992b) resuhs. 

Although it is obvious that the scatter in samples con- 
taining spiral galaxies is significantly larger than in sam- 
ples containing only E/SO galaxies (see also Visvanathan & 
Sandage 1977, Peletier & de Grijs 1998), we have shown 
that by avoiding the disturbing effects of the in-plane dust 
lane the scatter in the spiral galaxy NIR CM relation can 
be reduced significantly (e.g., Fig. |l|). 

Therefore, a NIR CM relation for spiral galaxies based 



5 IMPLICATIONS 

5.1 The shape of the CM relation 

Lasker (1970) was the first to notice that the CM relation 
seems to flatten towards the brighter galaxies in his sample, 
although in Fig. 14 of de Vaucouleurs (1961) one can already 
see this effect. Frogel et al. (1978) and TuUy et al. (1982) also 
note that the V-K colour of their brightest sample galaxies is 
relatively independent of luminosity compared to the fainter 
groups. 

With this relative flattening in mind, a linear flt to the 
CM relation may not be the best representation of the cor- 
relation. However, flattening effects are second-order effects; 
to first order linear fits should be sufficient for the determi- 
nation of distance moduli. 

To investigate a possible flattening trend in our data 
we computed the mean I-K colours in absolute magnitude 
bins of 0.5 mag, of which the result is shown as the solid 
line in Fig. H; the error bars indicate the dispersion in the 
data points in each bin. Although we are dealing with rel- 
atively few data points, a flattening is indeed appreciated 
at the bright end of the relation. This means that the old 
stellar populations in the brightest - or largest - disc galax- 
ies have a roughly constant intrinsic I-K colour. Thus, for 
spiral galaxies with Mk < —25.5 the CM relation cannot 
be used as a diagnostic tool for distance determinations. At 
the faint end too few data points are available to draw firm 
conclusions about the shape of the CM relation, but it seems 
likely that the correlation is maintained down to (field) spi- 
ral galaxies as faint as Mk ~ —20. 



5.2 Applications of a tight CM relation 

As we argued in Sect. H, a NIR CM relation for spiral galax- 
ies based on dust-free colours may in principle be useful as 
a diagnostic to estimate distances: it provides independent 
distance estimates based on observational parameters. 

For distance determinations to individual spiral galax- 
ies, we have to take into account the scatter in and the 
shallow slope of the CM relation, as well as the observa- 
tional errors, which will result in relative distances accurate 
to ~ 35% for the observations presented in this paper. This 
accuracy estimate is based on the observed dispersion (i.e., 
the combination of observational and intrinsic scatter) in the 
NIR CM relation in the absolute K-band magnitudes (~ 0.7 
mag; i.e. the dispersion along the horizontal axis in Fig. hla). 
For a "typical" sample galaxy, the observational dispersion 
in the absolute magnitude determinations (taking into ac- 
count the uncertainties in the apparent magnitude and the 
correction for the motion with respect to the Local Group) 
is ~ 0.45 mag, which means that the intrinsic scatter in 
the NIR CM relation among galaxies is ~ 0.5 mag for the 
absolute Tf-band magnitude determinations. Although it is 
likely that part of this intrinsic scatter is caused by the pe- 
culiar velocities of our sample galaxies with respect to the 
Hubble-flow movement of the centroid of the Local Group, 
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Figure 3. Composite CM relation, derived from samples (b) and (c) and supplemented with the observations presented in this paper. 
Our observations are shown as filled circles; the open squares correspond to sample (b), the open triangles to sample (c). The best-fitting 
composite CM relation is indicated by the dashed line; its parameters are listed in Table eI The NIR CM relation derived from the 
observations presented in this paper is indicated by the dotted line. The solid line was obtained by binning the data points in Mp^ bins 
of 0.5 mag; the error bars indicate the dispersion in the data points in each bin. 



small-number statistics prevent us from reaching conclusive 
evidence for this: peculiar motions are available for only 15 
of our 22 sample galaxies (Table m Mathewson et al. 1992) 
For instance, correcting the absolute magnitudes for the ef- 
fects of the peculiar velocities results in a dispersion in Mk 
of 0.994 mag, compared to 0.640 mag before correction of 
these 15 sample galaxies. The dispersion in {I — K) increases 
from 0.102 to 0.161 mag, based on bivariate fits. 

The observational uncertainties in the apparent Tf-band 
magnitudes discussed in this paper are largely due to the 
relatively shallow ("snap-shot") nature of the observations, 
as well as to the rapidly varying sky background in the K 
band (see de Grijs [1997], and de Grijs et al. [1997] for a 
discussion of the relevant photometric accuracy achieved 
in this project). This rapid variability requires that short 
(< 60s) object and sky observations, of equal integration 
times, are taken alternately, in order to correct for these 
background variability effects. Obviously, a greater number 
of such observations will smooth out any remaining effects 
due to the varying background on the one hand, and increase 
the signal-to-noise (S/N) ratio in the galactic outskirts on 
the other, thereby improving the photometric accuracy of 
the integrated ii'-band magnitudes significantly. In addition, 
to unambiguously establish a uniform spiral galaxy CM rela- 
tion, the distances to the sample galaxies need to be known 
to high accuracy. For our present sample, uncertainties in 
the absolute K-hand magnitudes can be reduced by the de- 



termination of distances that are more accurate than those 
used in this paper, which are based on recessional velocities. 

Therefore, by obtaining high-quality observations of 
and improved distances to moderately or highly inclined 
galaxies, the accuracy of this method can be enhanced by 
reducing the observational scatter in the apparent K-hand 
magnitudes to a lower limit of ~ 0.02 mag. The accuracy will 
thus be limited by the intrinsic dispersion in the NIR CM 
relation in the absolute A'-band magnitudes (~ 0.5 mag). Be- 
cause of this intrinsic scatter, the maximum accuracy of dis- 
tance determinations to spiral galaxies that can be reached 
using this relation will therefore be ~ 25%. Note that this ac- 
curacy approaches the 20% accuracy of Bower et al. (1992b) 
for early-type galaxies, for which the intrinsic scatter in the 
CM relation is significantly smaller. 

The intrinsic dispersion in absolute /("-band magnitudes 
in the CM relation for field spiral galaxies is most likely 
caused by the non-negligible effects of different degrees of 
active star formation, star formation histories, ages, metal- 
licities, and extinction, even at the 2 heights at which we de- 
termined our "dust-free" colours, and perhaps non-negligible 
vertical colour gradients in the discs of our sample galaxies 
(Peletier & de Grijs 1998; de Grijs & Peletier, in prepa- 
ration). In view of the expected variation of all these pa- 
rameters, it is even more surprising to find a very tightly 
constrained NIR CM relation for field spiral galaxies. 
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5.2.1 How Hubble Space Telescope observations can 
contribute to the calibration 

Since the galaxies need to be spatially resolved, high- 
resolution observations done with the Hubble Space Tele- 
scope (HST) can provide a powerful tool to minimize the 
observational scatter and extend the useful distance range. 
Such observations may therefore provide the means to cali- 
brate the relation. We will demonstrate this in the following 



• Model galaxy. We assume that the intrinsic major-to- 
minor axis diameter of a "typical" spiral galaxy disc is 9 : 1 
(e.g., Guthrie 1992), and the scale height ratio of the dust 
to the stellar disc is 1 : 2 (following the standard models in 
Huizinga 1994; see also Xilouris et al. [1999] for case studies), 
whereas the dust disc is embedded in the stellar disc accord- 
ing to the so-called "triplex model" (e.g., Disney, Davies & 
Phillipps 1989). Furthermore, either disc is assumed to be 
of exponential form, both radially and vertically. 

• Inclination restrictions. The requirement that we need 
to be able to control - to some extent - the effects of ex- 
tinction restricts the useful inclination range to inclinations 
between ~ 80° and 90°. The lower inclination limit is set 
by the fact that at i ~ 80° the dust lane will be observed 
at the outer rim of the projected galactic disc nearest to 
the observer. This implies that, from this inclination up- 
wards, the upper stellar layer at the far edge of the disc will 
not significantly be affected by the effects of extinction. For 
higher inclinations, an increasingly large dust-free area will 
be available for colour measurements, divided in two parts 
by the high-extinction in-plane dust lane. 

• Size restrictions. For an unambiguous determination of 
these dust-free colours, one will ideally need a dust-free area 
of at least three resolution elements containing galactic emis- 
sion brighter than the limiting magnitude or S/N ratio. With 
the superb angular resolution of the HST, we will be able to 
reduce the minimum projected sizes of the useful candidate 
galaxies considerably. In the NIR, the most suitable HST 
instrument is the Near-Infrared Camera and Multi-Object 
Spectrometer (NICMOS), which has a pixel size of 0."075 
(camera 2) and a FWHM for point sources of ~ 1.3 pixels 
(0."10) at l.l^m to 2.5 pixels (0."19) at 2.2/im (Krist & 
Hook 1997). The inclusion of a red optical passband (e.g., 
/ ) requires additional imaging with the Wide Field Plane- 
tary Camera 2 (WFPC2); the pixel scale of the Planetary 
Camera is 0."0455, with a point source FWHM of about 
2.5 pixels (~ 0."ll, depending on wavelength). Since the 
WFPC2 characteristics fall inside the parameter range for 
NICMOS camera 2, we will restrict ourselves to the possi- 
bilities provided by the use of HST NICMOS observations 
in the remainder of this exercise. 

Along the galaxies' minor axes, the useful area for colour 
determinations is, for all inclinations (80° < J < 90°) and as- 
suming a minimum dust-free extent of 3 resolution elements 
on the side nearest to the observer, ~ 0."29 and ~ 0."56 for 
the J and K bands, respectively. Under the additional as- 
sumptions outlined at the beginning of this paragraph, and 
by using Hubble's (1926) relation between a galaxy's inclina- 
tion and its axis ratio, the minimum projected size of a can- 
didate galaxy in the J band is l."6 x 8."5 to l."2 x 10. "5 for 
inclinations of 80° to 90° ; in the K band the different FWHM 



converts to projected sizes of 3."0 x 16."4 to 2."3 x 20."2 for 
the same inclination range. We are thus limited by the sizes 
required by the K-hand imaging. 

Because of the high resolution of the HST one will be 
able to resolve a galaxy of a given size (determined by the 
limiting magnitude or S/N ratio of the observations) out to 
larger distances than when observed with (seeing-limited) 
ground-based telescopes. To get an indication of the useful- 
ness of the NIR CM relation as a function of distance we 
obtained HST archive observations of NGC 891. This edge- 
on Sbc galaxy is comparable in size and appearance to the 
Galaxy. We used observations in the optical F606W pass- 
band (160s) and the NIR F160W filter (192s). 

In the (wide) V-band observations, the maximum useful 
dust-free region on the minor axis, where the galactic V- 
band surface brightness Ey < 25 mag arcsec"^, is (33. "0 ± 
0."5), which corresponds to (1.50 ±0.02) kpc at its distance 
of 9.5 Mpc (van der Kruit & Searle 1981, see also Keppel 
et al. 1991). The minimum dust-free region on a galaxy's 
minor axis containing three independent resolution elements 
is 0."56, which is required by the need for NICMOS X-band 
imaging, as explained above. 

From this argument it follows immediately that we can 
use galaxies of similar size as NGC 891 for the construction 
of a NIR CM relation up to distances of (560 ± 5) Mpc. 
For comparison, ground-based observations can be used for 
similar galaxies out to distances of only ~ 210 Mpc, un- 
der optimal seeing conditions (FWHM ~ 0."5). In Fig. 
we show the difference between the results obtained using 
HST observations and ground-based observations taken un- 
der such optimal seeing conditions for edge-on galaxies like 
those presented in this paper. To construct this figure, we 
have adopted the assumptions about the galactic structure 
outlined above, and a major axis diameter of D25 ~ 13. '5 
for NGC 891 (RC3), corresponding to a linear size of 37.3 
kpc. 

From a comparison between the optical and the NIR 
HST observations of NGC 891, with integration times of 
the same order, it is clear that the NIR observations do not 
reach the galactic outskirts at similar S/N ratios as in the 
optical. Therefore, if we base our conclusions on this par- 
ticular case, the successful use of the NIR CM relation for 
distance estimates depends on the faintest surface brightness 
levels that can be achieved in the NIR. However, this does 
not need to be the case for all galaxies, as it depends on pa- 
rameters, such as intrinsic colours, sky background (which is 
supposed to be very low for space-based NIR observations) , 
etc. 

In summary, high-resolution observations done with the 
HST can significantly increase the number of useful galax- 
ies and hence provide a powerful tool to calibrate the CM 
relation for spiral galaxies, reduce the observational scatter, 
and extend the useful distance range by more than a fac- 
tor of 2 compared to ground-based observations. The use- 
ful projected galaxy sizes are limited by the i<'-band resolu- 
tion, which translates to minimum sizes from 3."0 x 16. "4 to 
2."3 X 20. "2 for inclinations of 80° to 90°. Obviously, / (or 
J) and K-hand observations of galaxy fields from the HST 
data archive (e.g., the Hubble Deep Fields, among others) 
would be well suited to undertake such a programme. 



© 1999 RAS, MNRAS 000, 000-000 



10 R. de Grijs and R.F. Peletier 




20 30 

Diameter (kpc) 



-26 



-24 



-22 



-20 



-IE 



o J 

•K 



iSi 



iX'i 



'- ^.isi = -7.99 (±0.17) log (v^„j^,J - 2.69 (±0.36) . 
K: M^^^jf = -8.09 (±0.52) log (v^^^^^J - 4.69 (±1.09) 



100 



Figure 4. Comparison between the maximum distances out to 
which the NIR CM relation can be used successfully to obtain 
distance estimates for HST and ground-based observations under 
optimal seeing conditions. The power of space-based observations 
for distance determinations using this relation is clear. 



6 THE NIR TULLY-FISHER RELATION 

For the majority of our sample galaxies a homogeneous data 
set containing both the maximum (optical) rotational veloc- 
ities and total (apparent) /-band magnitudes was provided 
by Mathewson et al. (1992) and Mathewson & Ford (1996). 
In de Grijs (1997, 1998), we compared our /-band photom- 
etry to that of Mathewson et al. (1992) and Mathewson & 
Ford (1996). From the detailed comparison of our photome- 
try to theirs it was shown that we can accurately reproduce 
their results ({m/,our — m/, Mathewson) = —0.07 ± 0.13; de 
Grijs 1997, 1998). 

The / and /('-band TFRs derived from our absolute mag- 
nitudes, and supplemented with Mathewson's rotational ve- 
locity data, are shown in Fig. BJ We did not correct our 
integrated magnitudes for the effects of internal interstellar 
extinction; inclination corrections that aim to correct inter- 
nal extinction are generally not applicable to highly inclined 
galaxies. Moreover, the amount of extinction in a galaxy 
varies as a function of galaxy type, as well as among galax- 
ies of the same type (see, e.g., de Grijs et al. 1997, de Grijs 
1997, 1998, and references therein). It is therefore likely that 
the scatter in the / and K-hand TFRs derived from our ob- 
servations is at least partly caused by internal extinction. 
However, note that errors in the TF distances are generally 
dominated by uncertainties in the inclinations of the sample 
galaxies. 



6.1 The scatter in the TuUy-Fisher relation 

A good measure of the accuracy of the TFR can be ob- 
tained by assessing the scatter in the relation. In the / and 



Figure 5. /and if-band TFRs derived from our observations (no 
internal extinction correction was applied). 



K bands, the TFR r.m.s. scatter is 0.145 and 0.296 mag, 
respectively, based on bivariate fits (e.g., Giovanelli et al. 
1997, Verheijen 1998). The scatter in the /('-band TFR is 
significantly larger than that in the /band. This is probably 
caused by the fact that both the S/N ratio in our K'-hand 
images is significantly smaller and the effects of the varying 
sky background are much greater in K' than in the /-band 
observations, causing the larger uncertainties, and thus ob- 
servational scatter, in the K-hand magnitudes. We are con- 
fident that the scatter in the K band TFR due to internal 
extinction is very small in comparison, since the K-hand ver- 
tical surface brightness profiles exhibit hardly any signatures 
of either a regular dust lane or a patchy dust distribution, 
even in those galaxies showing pronounced dust lanes in the 
optical passbands (de Grijs et al. 1997). 

The r.m.s. scatter in our /-band TFR is small compared 
to that found in previous studies (for a discussion see Bern- 
stein et al. 1994, Giovanelli et al. 1997): only Bernstein et al. 
(1994) reported a significantly smaller scatter, of ct = 0.10 
mag, based on a sample of 22 galaxies. Giovanelli et al. 
(1997) argue that TF fits with scatter smaller than 0.25 mag 
are likely to be statistical accidents, which can occur when 
galaxy samples are small. They base this argument on the 
detailed study of the TFR for a large number of galaxies in 
24 clusters. Giovanelli et al. (1997) determined the scatter 
in the /-band TFR for the individual clusters; their obser- 
vational scatter varies between 0.12 and 0.36 mag (based on 
bivariate fitting routines) . With this statistically significant 
result in mind, we believe that the low scatter found in this 
paper is indeed a statistical accident. Additional support for 
this assumption is lent by the fact that the scatter in the CM 
relation increases (in either of the quantities) if we correct 
for the effects of peculiar motions, as shown in Sect. [3.4 

To our knowledge, studies of the TFR in the NIR K 
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(or K') band are scarcely available to date. From Fig. 10 
of TuUy & Verheijen (1997), we derive a scatter of 0.531 
mag in the ii''-band TFR of Ursa Major cluster galaxies, 
based on a bivariate fit; for the low and high surface bright- 
ness galaxies in this figure, we derive a scatter of 0.654 and 
0.363 mag, respectively. In addition, in a detailed study of a 
volume-limited complete sample of ~ 40 Ursa Major cluster 
galaxies Verheijen (1997) finds that the smallest observed 
TFR scatter in his galaxy sample is 0.29 mag in K' , consis- 
tent with no intrinsic scatter. 

Giovanelli et al. (1997) applied a type-dependent cor- 
rection for internal extinction to their observations, which 
varies between 0.50 and 1.00 mag in the /band, correspond- 
ing to 0.12 - 0.23 mag in the K band (assuming an Rv =3.1 
extinction law, e.g., Schlegel et al. 1998). If we apply a simi- 
lar extinction correction to our data, we find a scatter in the 
/-band TFR of 0.225 mag (as opposed to 0.145 mag for the 
uncorrected measurements), and a scatter of 0.316 mag in 
the K-hand TFR (compared to 0.296 mag before extinction 
correction), which is entirely consistent with Giovanelli et 
al.'s (1997) results. 

From their detailed analysis of the scatter in the TFR, 
Giovanelli et al. (1997) concluded that the average total scat- 
ter in their measurements of a statistically significant num- 
ber of galaxies in 24 clusters is ~ 0.35 mag. If we assume 
that this is a typical value for the scatter in any well-sampled 
TFR, the accuracy of the TF distance to a "typical" galaxy 
in our sample would amount to ~ 15%. 

When comparing our TF results with the accuracy of 
distance determinations using the NIR CM relation, we note 
that, due to the shallower slope of the CM relation, the ef- 
fects of uncertainties in the distance estimates are relatively 
more important for the cahbration of the spiral galaxy CM 
relation than for the TFR. On the other hand, one of the 
main advantages of the use of the NIR CM relation is that 
we only need photometric data to obtain distance estimates, 
whereas use of the TFR requires additional kinematic data. 



6.2 



Distance estimates: Tully-Fisher vs. 
Colour-Magnitude 



In Sect. 5.2.1 we estimated that the NIR CM relation can 



be used to obtain distances to galaxies similar to the Galaxy 
or NGC 891 up to (560 ± 5) Mpc, if HST-based photometry 
and resolution can be obtained. Following similar lines of 
argument, in this section we will estimate the correspond- 
ing maximum TF distance to a Galaxy-type system that 
can still be determined unambiguously. To do so, we will 
restrict ourselves to ground-based rotation curve measure- 
ments, since they are most readily available in the literature, 
and most easily obtained. 

Based on a statistical analysis of a spiral galaxy sample 
designed for TF applications, Courteau (1997) has shown 
that the best measure of TF velocity is given at the loca- 
tion of peak rotational velocity of a pure exponential disc. 
Alternatively, one can use the 20% width of the Hi veloc- 
ity profile. Although this latter method does not make any 
a priori assumptions about the luminosity profile or shape 
of the rotation curve, its accuracy is not optimal (Courteau 
1997). 

If we restrict ourselves to optical observations, one needs 
at least 5 resolution elements with reliable S/N ratios to ob- 



tain an accurate measure of the rotation curve. Under opti- 
mum seeing conditions (FWHM ~ 0."5), this corresponds to 
a minimum angular diameter of a galactic disc of ~ 2. "5. A 
Galaxy-type sy stem like NGC 891, with a linear diameter of 
37.3 kpc (Sect. 5.2.1), would therefore be useful for distance 
estimates up to ~ 3.1 Gpc (~ 5.5 times as far as the maxi- 
mum CM distances). Obviously, the use of Hi observations 
will increase this distance estimate significantly. 

The main disadvantage of the use of the TFR versus 
the CM relation for distance estimates to relatively nearby, 
highly inclined galaxies is that the conversion of apparent 
to face-on corrected magnitudes is controversial and model 
dependent. 



7 SUMMARY AND CONCLUSIONS 

In this paper we have looked at two methods of distance de- 
termination to spiral galaxies in the field using optical/NIR 
observations, the (I-K) vs. Mk colour - absolute magnitude 
relation and the /and /f-band Tully-Fisher relation. We have 
used "dust-free" colours to achieve greater accuracy for dis- 
tance determinations using the CM relation compared to the 
integrated galaxy colours that are generally used for this ex- 
ercise. Our main conclusions are the following: 

• Our data, supplemented with observations taken from 
the literature form a well-constrained composite spiral 
galaxy CM relation; it appears that the NIR CM relation 
for the old-disc population of spiral galaxies in the field is 
universal. 

• By avoiding the disturbing effects of the in-plane dust 
lane the observational scatter in the CM relation can be re- 
duced significantly. Therefore, the NIR CM relation for field 
spiral galaxies, based on dust-free colours, may in principle 
be useful as a diagnostic tool to estimate distances with an 
accuracy of ~ 25%. This accuracy is limited by the intrinsic 
dispersion of the NIR CM relation in the K-h&nd absolute 
magnitudes (~ 0.5 mag). 

• High-resolution observations done with the Hubble 
Space Telescope can provide a powerful tool to reduce the 
observational scatter, calibrate the relation, and extend the 
useful distance range (by more than a factor of 2 compared 
to ground-based observations). The useful projected galaxy 
sizes are limited by the NICMOS /f-band resolution, which 
translates to minimum sizes from 3."0 x 16. "4 to 2."3 x 20. "2 
(at the limiting magnitude or S/N ratio) for inclinations 
80° < i < 90°. 

• We found a scatter in the / and /f-band TFRs of 0.145 
and 0.296 mag, respectively. Although, to date, these values 
are among the lowest found in these passbands, they are 
likely statistical accidents, caused by our relatively small 
sample size (see Giovanelli et al. 1997). 

• Due to the shallower slope of the CM relation, the scat- 
ter is relatively more important for the determination of 
distances via this method than for distance determinations 
using the TFR: typically, TF distances can be determined 
with an accuracy of ~ 15% (and to significantly greater dis- 
tances), as opposed to the ~ 25% accuracy achieved using 
the NIR CM relation. However, one of the main advantages 
of the use of the NIR CM relation is that we only need pho- 
tometric data to obtain distance estimates, whereas use of 
the TFR requires additional kinematic data. 
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• Although a linear fit is a good first-order approximation 
to the composite NIR CM relation, the old-disc populations 
of the brightest disc galaxies {A'Ik < —25.5) have a roughly 
constant intrinsic I-K colour. Our observations at the faint 
end of the CM relation are consistent with a linear correla- 
tion down to galaxies as faint as Mk ~ —20. 
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